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ABSTRACT 

Annua] extremes of total column amount of ozone 
(in the period 1964-1988) from a network of29 Dobson 
stations have been examined using the extreme value 
analysis. 

The extremes have been calculated as the highest 
deviation of daily mean totalozone from its long-term 
monthly mean, normalized by the monthly standard 
deviations. The extremes have been selected from the 
direct-Sun total ozone observations only. The 
extremes resulting from abrupt changes in ozone 
(day to day changes greater than 20%) have not been 
considered. 

The ordered extremes (maxima- in ascending way, 
minima - in descending way) have been fitted toone 
of threeforms of the F i s h e r- T i pp e t extreme value 
distribution by the nonlinear least square method 
(Levenberg-Marguard method). 

We have found that the ordered extremes from a 
majority of Dobson stations lie close to Fisher-Tippet 
type III. 

The extreme value analysis of the composite annual 
extremes (combined from averages of the annua] 
extremes selected at individual stations) has shown 
that the composite maxima are fitted by the Fisher- 
Tippet type III and the composite minima by the 
Fisher-Tippet ty pe I. The difference between the 
Fisher-Tippet types of the composite extremes seems 
to be related to the ozone downward trend. 

Extreme value prognoses fort he period 1964-2014 
(derived from the data taken at: allanalyzed stations, 
the North American, and the European stations) have 
revealed that th e progn ost ic extremes are close to 
the largest annual extremes in the period 1964-1988 
and there are only sm a 1 1 reg i on a I differences in the 
prognoses. 

1. INTRODUCTION 

Analyses of both ground-based [IOTP, 1990] and 
satellite ozone data [Stolarski et al., 1991 ] indicate 
statistically significant total ozone decline. 

The potential effects of an increase in solar 
ultraviolet radiation on the biosphere have 
stimulated m uch of the concern about ozone 
depletion. Then, the estimation of future ozone 
changes is of special interest for human life. 

In this paper we examine possibilities to implement 
the theory of extreme value toinfer future total 
ozone extremes. Many extreme value analyses of 


hydrological and meteorological variables have 
been undertaken in the past (for review, see 
Tabony, 1983). The method provides estimation of 
the highest and the lowest values which 
meteorological (orother) variable is likely to attain 
in a given number of years. The extreme value 
prognoses are based on the assumption that all 
factors determining changes in a prognostic 
variable have already appeared during theperjod of 
observations. 

2. DATA 

Total ozone extremes inthe period 1964-1988 have 
been selected from a network of 29 Dobson stations. 
The names of the stations, index in the WorldOzone 
Data Center, locations and the data periods arc listed 
as follows: 


1 Lerwick, 

U.K 

43, 60°N, 1° W ; 

1/64- 1 1/88 

2 St. Peters burg Russia 

42, 60°N 30°E; 

8/68-12/88 

3 Churchill, 

Canada 

77, 59°N, 94°W: 

l 65-12/88 

4 Edmonton, 

Canada 

21, 54°N, 1 1 4°W; 

1/64-12/88 

5.Belsk, 

Poland, 

68, 52°N, 2 1° E ; 

1/64- 1 2/88 

6 Bracknell, 

U.K. 

102, 5 1°N, 1<HV; 

1/69-12/88 

7 Uccle, 

Belgium 

53, 5 1°N, 4°E; 

2/71-12/88 

8 Hradec K., 

Czech. 

96, 5 0° N , 1 6° E ; 

1/64-12/88 

9 Hohenpciss 

, Germany 

99 48°N, 1 1°E; 

1/67-12/88 

10 Caribou, 

U.S.A. 

20, 4 7° N , 68° W; 

1, 64- 1 2,88 

11 Arosa, Switzerland 

35, 4 7°N , 10°E; 

1/64-12,88 

12 Bismarck, 

U.S.A 

19, 47° N , 1 0 1° W ; 

1/64- 1 2/ 88 

1 3 Toronto, 

Canada 

65, 44° N, 7 9° W ; 

1/ 64-12,88 

14 Sapporo, 

Japan 

12, 43° N, 1 4 1°E; 

1 64-12,88 

15 Rome, 

Italy 

55, 42°N, 1 2° E ; 

1/64- 1 2,88 

1 6 Boulder, 

U.S.A. 

67 40°N , 1 05° W ; 

1/64-12,88 

17 Cagliari, 

Italy 

38, 3 9° N , 9°E; 

1,64-12,88 

18 Wallops Is., 

U . S : A . 

107, 3 8° N , 76° W; 

l, 70-12 88 

19 Nashville, 

U.S.A. 

106, 36° N, 8 7° W ; 

1/64-12 88 

20 Tateno, 

Japan 

14, 3 6° N , 1 4 0° E ; 

1.64-12 88 

21 Srinagar, 

India 

13, 3 4° N , 75°E; 

2/64-12 88 

22 Kagoshima, 

Japan 

7, 32°N, 1 3 1°E ; 

1, 64-12 88 

23 Quetta, 

Pakistan 

11, 30°N, 67°E ; 

8 69- 1 2. 88 

24 Cairo, 

Egypt 

152, 30°N, 3 l 0 E : 

t 1, 74-12/8 8 

25 Naha, 

Japan 

190, 26°N , 1 28°E 

4, 74- 1 2, 88 

26 Mauna Loa, 

U.S.A. 

31, 19°N, 1 5 6° E ; 

1 64-12 88 

2 7 Huancavo, 

Peru 

110, 1 2° S , 75° W ; 

2. 64- 1 2, 88 
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28 Samoa. U.S.A. 191, 14°S, 171°W; 1/76-12/88 

29 Aspendale , Australia 26, 38°S, 145°E; 1/64-12/82 

For each station, the annual extremes in total column 
amount of ozone have been calculated fromthe daily 
averages of total ozone published in Ozone Data forthe 
World (ODW) journals. However, International Trend 
Panel-1988 [IOTP, 1990] have found strong 
incongruities in the published data of many ground- 
based stations due to instrument calibrations. 

A provisionally revised data set (1957-1986) was 
prepared by the Panel, with t h e c o r r e c t i o n s being 
applied to the monthly averages of total ozone values. 

Making use of thePancl's re-evaluated total ozone 
monthly means we convert all daily averages of total 
ozone values (calculated before Jan. 1987) to the 
provisionally corrected ones in the following way; 

. . _ M l(j,k,l) 

03 n e w ( t ,j , k 1 1 ) - old ( 1 » J * k * * ^ M 2 (j , k , l ) 

where: O 3 n € w (ij,k f lh 0 3 0 / d (iJ,k,D ‘ daily mean 
of total ozone o n * - 1 h d a y , j - 1 h month, k - 1 h year and 
at / -th Dobson station, provisionally revised and 
unrevised, respectively, M I ( j , k , I )- re vised monthly 
mean oftotalozone at / station in 7-th month and /c - 
th year (data from IOTP, 1990), M 2 (j , k , l ) - monthly 
mean of total ozone derived from ODW data. 

Provisionally corrected daily a v e r a g e s of total ozone 
values are transformed to the departures from the 
long-term monthly means ( l 964- 1988) expressed in 
percents of the long-term monthly standard 
d e v i a t i o n . F i 11 a 1 1 y , the annual extremes have been 
selected from the dimensionless daily averages o f t o t a 1 
ozone values. 

ODW data from only a few stations (Belsk, Sapporo, 
Kagoshima, Tsukuba, Naha) could be accepted without 
any adjustment of the data, because original reported 
results of total ozone measurements were re-evaluated 
(and published in ODW) on daily basis taking into 
a c c 0 u n t a 1 1 instrument's calibrations. Our correction 
method provides only a crude estimation of daily 
mean total ozone. 

ODW data contain results from the various types oi 
o/onc observations among these the dircct-Sun oz o n e 
observations provide the highest accuracy. Taking 
this into consideration, the e x t r e m e s h a v e been 
selected only fromthe direct-Sun observations. The 
e x t r e m e s , w h i c h were formed as a result of abrupt 
changes in t h c ozone pattern (day-to-day change in 
1 0 1 a 1 ozone more than 2 0 T; ) , h a v e not b e e it c o 11 s i d e r e d 
Therefore, the extremes i n t o t a 1 ozone, examined in 
this paper, are caused b y p r o c e s s e s w i t h a time scale 
larger t h a n a few days. Moreover, in the case who 11 
the day-to-day change in total ozone i s v e r v large 
there is a larger likelihood that this change is a result 
of erroneous observations. 

3 . THE EXTREME VALUE ANALYSIS 

The t h e o r y of extreme value was first developed b y 
Fisher and Tippet, [ l 928 ] . They showed that extremes 
taken from a sample containing N observations 
converge asymptotically towards one o foul y t li r e e 
forms as N increases. These are called Fisher-Tippet 
Type 1, Type II and Type 111. A general solution of the 
extreme value problem was obtained b y J e n k i 11 s o n 
[ 1 955 ) in the form ; 


X = X 0 + a 


1 -exp(-fc Y) 
k 


[ 1 ] 


where: X- extremes (ordered in ascending way), 

Y - reduced variate calculated from the equation: 

F(*) = exp( -exp(-Y) )> 

F(AQ- empirical cumulative probability function, X G 
a , andk constants to be calculated. 

Fisher-Tippet types I, II, and III are characterized 
by their different shapes when plotted on a graph X 
against Y (see Fig. 1). Type Icorrespondsto* = 0 and 
forms a straight line. Type II (bounded below but 
not above) a n d T y p e III (bounded above but not 
below) correspond to k < 0 and k > 0, respectively. 



Fig.l. Fisher-Tippet distribution types I, //, /// 

4. RESULTS 

I n order to use extreme value method, the 
empirical cumulative distribution function of 
annual extremes (ECDFAE) has to be constructed. The 
annual extremes are ordered fromthe lowest A] (for 
minima we use absolute values) to the h i g h e s t A m , 
where m is the number o t annual extremes. ECDFAE 
could be calculated b y various formulas, wehave 
chosen Beard's formula [Beard, 1943]; 

F (A/ ) = (1 - 0 . 3 1 ) -■ ( m + 0 . 3 8 ) 


F i n a 1 1 y , o n e of the three Fisher-Tippett types has to 
be fitted to a set of (A',\ V/) points, 

Yj= - 1 n ( - 1 n ( F (A/ ) ) ). 

In practice, after plotting the s e t o I (A 7, Y / ) points 
on agraph A a g a i 11 s t > , the F i s h e r - T i p p e t type of 
ECDFAE is d e t c r mi 11 e d " h y e y e " method and t h e 
parameters describing precisely the pattern o f 
ECDFAE are calculated assuming that the t y p e o f 
ECDFAE is known. However, in our paper we use the 
nonlinear least square method (Lev e 11 her g- 
Marguardt method) to find the model's constants; 

Xo . a 1 and fc . 

Then, both the values of A* and the standard error of 
k let us determine the type of ECDFAE. 

In Fig. 2 , the estimates of A from individual stations 
extremes data are presented. It is seen that Fisher- 
Tippet type 111 dominates. 
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Fig. 2. Values o f k for annual maxima (Fig. 2 A) and 
for annual minima (Fig. 2 B ) derived from the annual 
extremes taken at the selected Dobson stations. 

Number above each point corresponds to the ODW 
station index. Vertical bars represent standard errors 

When, the extremes lie close to the physically imposed 
bound, a type 111 distribution of the extremes maybe 
expected, but when the observed extremes fall far 
from their limit, they may appear to befitted by types 
I or II. However, type III and II may appear 
erroneously if the extremes are selected from small 
population (too small number of direct-Sun 
measurements per year) or the data contain samples 
from \ arious populations (extremes are caused by 
more than one mechanism ). Type Mean be form 
misleadingly if the data are contaminated by only a 
few v ery large extremes coming from erroneous 
observations. 

Then, t o determine the type of extremes from 
individual station data correctly we need to a n a 1 v z c a 
cause of the extremes appearance. 

Examination of the composite extremes, combined 
from the annual extremes taken at each station, is an 
attempt to eliminate the observation errors. In Fig. 3 
we present result of extreme value analysis applied to 
the composite annual total ozone extremes. 

Maxima are well fitted by Fisher-Tippet type III 
distribution, A: = 0.4 6 + -0.06. The minima (absolute 

values) lie close to Fisher-Tippet type I distribution, k 
= -0.02 + .< - 0.06. These findings imply that the future 
maxima with the values much larger than the largest 
maximum registered up to 1988 will n o t b e expected 
(the highest maximum in the period 1 9 64-1988 lies 
close to upper limit imposed by Fisher-Tippet III type) 
.while more negative minima will be recorded (total 
ozone minima seem to be unbounded). 

The difference between the Fisher-Tippet types 
obtained for the annual maxima and minima seems to 


be related to decreasing tendency in atmospheric 
ozone. The ozone extremes have been calculated as 
the highest annual departures from the long-term 
(1 964- 1 988) m e an s . In the future, larger maxima 
(deeper minima) than the ones recorded up to 1988 
will beunlikely (likely) because the negative trendy 
in ozone attenuate the maxima (amplify (he minima) 
relative to the ozone mean level in the period 19 64- 
1988. 



Fig. 3. Empirical cumulative probability func t ion of 
annual extremes (for minima we use absolute 
values) combined from the annual total ozone 
extremes taken at 29 Dobson stations . 

The extreme value method provides the possibility 
to estimate the highest and the lowest v alues, which 
any meteorological (or other) variable is likely to 
attain in a given number of years. The estimation is 
obtained from the extrapolation of F i s h e r - T i p p e t 
distribution at a reduced variate value, which 
corresponds with a chosen period of prognosis. 

Tab.l shows results of extreme value prognosis in 
the period 1964-2014 for all Dobson stations listed in 
section 2 and for 10 European and 8 North American 
stations. We examine both provisionally revised data 
and non corrected ODW data. 

Table 1 . Extreme value prognoses (expressed in the 
standard deviations units ) in the period 1964-2014 
for 29 Dobson stations and for the European and the 
North American stations. Standard errors of the 
estimates in parentheses. 


Class 


REVISED 

MIN MAX 


DAT A 

NON-REVISEI) 

MIN MAX 


All stat. -2.79(0.02) 3.53(0.04) -2.68(0.03) 3.55(0.04) 

Europe -2.80(0.03) 3.91(0.05) -2.81(0.04) 3.90(0.05) 


N.Amer. -2.82(0. 04) 3.47(0.04) -2.79(0.03) 3.53(0.05) 


Comparing the extreme v alue prognoses for the 
anal v zed groups of ozone stations we have found 
that; provisionally revised data and non corrected 
ODW data give similar result, there are small 
regional differences in the prognoses. Prognoses 
for minima (absolute values) are lower than for 
maxima (this reflects positive asymmetry of the 
parent distribution). 
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In the period 1964-1988, the largest maximum (about 
3.5 standard d e v i a t i o n ) a n d minimum (about -2.6 
standard deviation), calculated for the group of all the 
analyzed stations (see Fig. 3), d i ff er only slightly from 
the prognostic extremes for the period 1964-2014. 

We estimate that 1 standard deviation for the group 
of all the analyzed stati ons eq ua Is about 10 Dobsons 
Then, the predicted largest minimum (for the period 
1964-2014) is only a few Dobsons lower thanthat 
observed in the period 1964-1988. Therefore, a 
downward trend in total ozone (about 1.5% per decade 
for the region 3CPN-60°N in the period 1970-1988, 
Krzyscin, 1992) will influence only a little the future 
largest extremes. 

Almost all studies of the long-term variations in 
atmosphere ozone (especially the trend models) were 
limited to analyses based on the monthly meansor the 
annual means of total ozone. Our paper reveals that 
the extreme value analysis provides a useful tool for 
the detection long-term changes in atmospheric 
ozone. The reliably of climatological statements based 
on the analysis of ozone extremes is clearly de pe n de nt 
upon quality of the data. Then, elimination of 
spurious extremes from the data seems to be the most 
important part of the analysis. 
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